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Abstract. 

The coupled evolution of the spin frequency, core temperature, and r-mode amplitude of an accreting neutron 
star is calculated. We focus on those conditions that can produce persistent gravitational radiation from the r- 
mode. During X-ray quiescent phases of transient LMXBs, one may be able to identify the constant contribution 
of the gravitational wave emission to the spindown rate. Another signature is the r-mode contribution to the 
heating. 
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There is now strong evidence [4] that the fre- 
quency f — 1/P of burst oscillations in accreting 
low-mass X-ray binaries (LMXBs) asymptotes to 
the spin frequency of the neutron star, as sepa- 
rately measured in the coherent millisecond X-ray 
pulsar LMXBs. Moreover, the burst oscillations are 
often also sufficiently coherent to allow a determi- 
nation of the spindown (or spinup) rate P~ 1 dP/dt. 
The distribution of spin periods indicates that grav- 
itational radiation may be limiting the accretion- 
driven spinup of these neutron stars [1, 17, 4]. Here 
we investigate the possibility that steady spindown 
produced by the emission of gravitational radiation 
via an excited r-mode might be detectable during 
the quiescent phases of these X-ray transients. 



Governing relations 

Details of the physical conditions and evolution 
of neutron stars that have been spun up by accre- 
tion to a state in which gravitational radiation can 
be emitted from an r-mode, and references to its 
earlier development, are provided by Wagoner [18]. 
The dominant (I — m = 2) r-mode (velocity pertur- 
bations 5v driven by the Coriolis force) is charac- 
terized by an amplitude a ~ Sv(r = R)/CIR, where 
CI = 2tt/P and R are the angular velocity and ra- 
dius of the neutron star. With J*(M, CI) the angular 
momentum of the equilibrium neutron star, its to- 
tal angular momentum is J = J* + (1 — Kj) J c , where 
the canonical angular momentum is J c = —K c a J*. 
(These constants Kq ~ 1.) 



The angular momentum perturbation obeys the 
relation 

dJ c /dt = 2J c [{F g (M,CTj-F v (M,Q,T)] . (1) 

The gravitational radiation growth rate is F g = 
T^, 1 (ri/n c ) 6 , with T gr ~ 3 s for the neutron star 
model adopted and Cl c = (nG(p)) 1 / 2 . In what fol- 
lows, we shall neglect the dependence on the slowly 
increasing mass M. 

The viscous damping rate is 

F v S F bl {S n ,S s ,B, CI, T) + F hb {T h ,Cl, T) , 

where T is the core temperature, and B is the 
magnetic field strength in the core-crust boundary 
layer. We have modified the viscous and magnetic 
boundary layer damping rate (Ff,i) of Kinney and 
Mendell [9], with the splippage factors [11] S for 
the normal and superfluid components. The hyperon 
bulk viscosity damping rate (-F^b) due to n + n ^ 
n + A, etc. employs results of Lindblom and 

Owen [12] and Haensel, Levenfish, and Yakovlev [7]. 
(But very recently van Dalen and Dieperink [16] 
obtained a smaller rate.) 

Conservation of total angular momentum then 
requires that 

dJ/dt = 2J c F g + Ja M(t) - T B (M) . 

The first (negative) term is due to gravitational 
radiation, the second to accretion [at the rate M(t)], 
and the third is the magnetic torque (e.g. [14]) of the 
disk on the star. 



From these relations, we can obtain the spin evo- 
lution equation: 

(i)f = -m 3 F 9 + (l-K 3 )F v ]K c ai 

+ fc^M(i)-^, (2) 

where J*(M,Q) = <9J*/<9fi and j*(M,Q) = dJ«/dM. 

Finally, thermal energy conservation for the star 
gives 

— cdK = C(T)—^2E C F V (Q,T) 

+ K n {M) n c 2 -L v {T), (3) 

where the r-mode viscous heating is propor- 
tional to the rotating frame canonical energy 
E c = (K c /3)Ji,ila 2 . We have assumed that the 
thermal conductivity timescales are short enough 
to allow the use of a single spatially-averaged core 
temperature T. The inner crust nuclear heating is 
driven by an accretion rate (M)„, averaged over the 
nuclear time scale r„ ~ 1 year [20]. The constant 
K n mlx 1CT 3 . 

Comparison of observations of thermal emission 
from isolated neutron stars with computed cooling 
histories has led Kaminkcr, Yakovlev, and Gnedin [8] 
to propose that at the core temperatures of interest 
here (T~3x 10 8 K) the core (triplet) neutrons are 
normal, while the core (singlet) protons and inner 
crust (singlet) neutrons are superfluid. The neutrino 
lumionosity that we employ [18] is then due to di- 
rect Urea, modified Urea, electron-ion and neutron- 
neutron neutrino bremsstrahlung, and Cooper pair- 
ing of (inner crust) neutrons. 

For typical values of ft, T, and (M) < Msdd ~ 
3 x 10 _8 M Q /yr (with j a ~ y/GMR), we obtain three 
key time scales (gravitational radiation, cooling, and 
accretion): r g = 1/F g ~ 10 3 s , 

T ^T7^~ 103yr < r a .-J->10V. 

M J ) Ja(M) a 

These time scales are seen to be relaxation times 
of our evolution equations (1), (3), and (2), respec- 
tively. Here the accretion rate is averaged over the 
time scale r . 



Evolution 

We are mainly interested in the evolution of neu- 
tron stars after they have been spun up to the point 
where equation (1) vanishes: F g (Q in ) — F v (Q in ,T in ). 



This equality defines our initial state, where the per- 
turbation can begin to grow. The initial temperature 
Ti n is determined by the vanishing of equation (3), 
with the nuclear heating balanced by the neutrino 
emission. 

In contrast to the initial state, the equilibrium 
state is defined by the vanishing of the evolution 
equation (2), in addition to the evolution equations 
(1) and (3). The equilibrium amplitude is then given 
by a eq = [r g /(2K c r a )} 1 / 2 ~ (1(T 7 - 10" 5 ), assuming 
that the magnetic torque is negligible. 

In Fig. 1, we show the critical curve defined by 
F g (Q) = F V (Q,T), for a somewhat optimistic choice 
of damping rate parameters. On it are indicated 
the initial and equilibrium states for two choices 
of (time averaged) mass accretion rate. The linear 
perturbation analysis of Wagoner, Hennawi, and Liu 
[19] shows that stability of the equilibrium state 
requires that the slope dft/dT > 0. 
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FIGURE 1. The relation between dimensionless an- 
gular velocity fi/fi c (w 2/3 at the shedding limit) and 
temperature T% on the critical curve defined by F g (Q) = 
F V (Q,T). The model chosen has a normal fluid hyperon 
damping rate 10 times greater than that of Lindblom 
and Owen [12], a hyperon superfluid transition temper- 
ature T h = 1 x 10 9 K, core— crust slippage factors S n — 
S s = 0.2, and boundary-layer magnetic field B < 10 9 G. 
Also shown are the initial and equilibrium states for (a) 
(M) a = Mfidd/300: T 8 (m) = 1.84, T 8 (eq) = 2.44 and for 
(b) (M) a = M Edd /3: T 8 (m) = 3.32, T s (eq) = 3.73. 

It can be shown that the evolution from the initial 
state is also controlled by the sign of the slope dQ jdT 
of the critical curve: 

• If the slope is negative [as in Fig. 1, case (a)], 
there will be a thermal runaway (dT/dt > 0, dfl/dt w 

1 /2 

0) with a growth rate K r ' a eq /r g <~ 1/yr that is 
of the same magnitude as found by Levin [10]. 
Here K r ~ 10 5 is the ratio of rotational to thermal 
energy in the star. The initial phase of the runaway 
occurs at roughly constant angular velocity. The 
later evolution is uncertain, but the neutron star 
may again reach the critical curve, where its slope 



could be positive. 

• If the slope is close to zero, there will initially be 
overstable oscillations of the type found by Wagoner, 
Hennawi, and Liu [19]. 

• If the slope is positive [as in Fig. 1, case (b)], the os- 
cillations of the growing amplitude are damped out 
on a timescale ~ r c , after which it slowly increases 
toward its equilibrium value, as shown in Fig. 2. The 
evolution of ft and T is similar, with the time re- 
quired to reach the equilibrium state given by 

At«[(fte,-n <n )/fi]T . 

If the magnetic torque term in equation (2) were 
important, r a could be appropriately modified. Af- 
ter accretion ceases, the star spins down along the 
critical curve. 




Log 10 (t/t„) 

FIGURE 2. The early evolution of the r-mode ampli- 
tude a, for M = MEdd/3. The initial amplitude (during 
spinup to the critical curve) was chosen to be Qo = 1CP 12 , 
and t* = 1 yr. When equilibrium is reached, a eq — 3.60 x 
1(T 6 . 

The ratio of the spindown rate to spinup rate is 



P- 1 dPM(M<(M) a ) 
P- 1 dP/dt(M=(M) a ) 
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We also find that on the critical curve, the ratio of r- 
mode heating to accretion-induced nuclear heating 
in equation (3) is 

8 / \ 2 
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In all cases, if the neutron star reaches a positive 
slope section of the critical curve in any way, equa- 
tions (1) and (3) quickly relax to equilibrium. So 
further evolution occurs at a rate governed by equa- 
tion (2) (whose first term is now —2K c F g a 2 ), with 
steady gravitational radiation generated by the r- 
mode amplitude given by the vanishing of equation 
(3) (with F g = F v ): 

a 2 = (L„ -K n (M) n c 2 )/[(2K c /3)J^F g ] . (6) 



This remains true as long as <Kl/dT > 0. After 
accretion ceases, the star will spin down until this 
slope vanishes [15]. Then the star will cool further at 
constant ft, no longer emitting gravitational waves 
(a negligible). It then can become a millisecond 
radio pulsar [15]. 

After T increases slightly above T in , the neutrino 
luminosity dominates the nuclear heating (for con- 
stant (M)n) in equation (6). We plot the resulting 
amplitude at any point on the positive slope section 
of the critical curve in Fig. 3. Using this amplitude, 
we show in Fig. 4 the spindown rate when the ac- 
cretion and magnetic torques are smaller than that 
due to gravitational radiation in equation (2). 
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FIGURE 3. The r-mode amplitude a at any point 
on the stable section of the critical curve where L v 3> 
K n {M) n c 2 . 
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FIGURE 4. The spindown rate produced by gravita- 
tional radiation when the neutron star is on the stable 
section of the critical curve. 

Observations of transient LMXBs during their 
quiescent states (when M(t) <C (M)) provide the op- 
portunity to detect two signatures that the neutron 
star is emitting gravitational radiation at the rate 
predicted if it has evolved to the stable portion of 
the critical curve. In order that such a positive slope 
portion exist, it appears that the neutron star core 
must contain hyperons with a moderate superfluid 



transition temperature (<~ 10 9 K) and neutrons with 
a low supcrfluid transition temperature (~ 10 8 K) 
[18]. 

The first signature was pointed out by Brown and 
Ushomirsky [3] , who showed that observations of the 
X-ray luminosity during quiescence can place inter- 
esting limits on the amount of r-mode heating. We 
have also estimated the photon luminosity and ob- 
tained results similar to theirs for a 'normal core'. 
They assumed that the star was in the equilibrium 
state. Note however from equation (5) that the mode 
heating will often be less than they obtained dur- 
ing spinup to or spindown from this state. Within 
the uncertainties of the mass and neutrino luminosi- 
ties consistent with observations of isolated neutron 
stars [21], comparison with observations of the four 
quiescent sources with large spin frequencies allows 
the presence of the predicted level of r-mode heat- 
ing. However, this extra source of heating is certainly 
not yet required [21]. 

Here we focus on the second signature. When the 
accretion torque becomes negligible at sufficiently 
small values of M(t), the constant gravitational wave 
torque provides a fixed lower limit to the spindown 
rate, shown in Fig. 4. If the contribution of magnetic 
torque can be identified from its dependence on M(t) 
and removed, a constant residual will provide evi- 
dence that the neutron star is emitting gravitational 
radiation at a significant rate. The transient XTE 
J0929-314 was observed to have an average spindown 
rate P~ 1 dP/dt w 1.6 x 1CP 8 per year [6], comparable 
to that in Fig. 4 for a core temperature T w 3.5 x 10 8 
K. The transient SAX J1808.4-3658 was observed 
to spin down at about the same rate after an out- 
burst in 2002, but spun up by a larger amount during 
its 1998 outburst [13]. As more data accumulate, it 
may be possible to extract more information about 
the contributions to spindown in the candidate (high 
spin frequency) neutron star transients. 

If Sco X-l has been spun up by accretion to 
a stable equilibrium state (in which gravitational- 
wave flux is proportional to X-ray flux) [17], it 
should be detectable by Advanced LIGO. When 
signal recycling (narrow-banding) is employed, a 
few additional LMXB's may also be detectable [2, 
5]. A major uncertainty is the actual value of the 
(very) long-term average (M) a , as inferred from 
the X-ray flux. If the neutron star is not in the 
equilibrium state, but is on the stable portion of 
the critical curve, the gravitational-wave amplitude 
(metric perturbation) will scale as h <x fl 3 a. For the 
dominant r-mode, the gravitational wave frequency 
is/ su ,^4/(3P) = 2Q/(37r). 
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